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I ■ Abstract. Fractional Levy motion (fLm) is the natural generalization of fractional 

^ I Brownian motion in the context of self-similar stochastic processes and stable 

^ ' probability distributions. In this paper we give an explicit derivation of the propagator 

of fLm by using path integral methods. The propagators of Brownian motion and 
fractional Brownian motion are recovered as particular cases. The fractional diffusion 
equation corresponding to fLm is also obtained. 



o 
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^ ' 1. Introduction 

> 

OO . 

. It is widely known that if ^2{t) is a Gaussian, uncorrelated noise (i.e. white noise) the 
Langevin (stochastic) equation 



O' x{t) = xo+ I Ut')dt', (1) 

00 ! Jo 

^ ' describes ordinary Brownian motion. One of the properties of Brownian motion is that 
the average squared displacement grows linearly with time, {{x{t) — xq)^) oc t. However, 
^ I many transport processes in physical, biological and social systems exhibit anomalous 
; diffusion [H [2]. That is, {{x{t) - xof) oc t^^ , with H ^ 1/2, where H is called the 
Hurst exponent [3]. In the past, several authors have attempted to generalize Eq. ([T]) 
in order to accommodate these anomalous processes. The anomalous behaviour may 
be associated to the existence of spatiotemporal correlations that produce correlated 
increments dx{t) := x{t + dt) — x{t). For that reason, the first proposed generalization 
had the form 

x{t) = Xo+ ^^hIi/2) ly ~ ^')''"'^'^2(0dt', (2) 

known as fractional Brownian motion (fBm) [4J, which has been extensively studied 
and applied [HI [6l [71 [8] . In Eq. ([2]) the noise ^2(^) remains Gaussian and uncorrelated 
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but, thanks to the convolution with the power-law kernel, the increments dx{t) become 
correlated in a way that yields the desired scaling of the mean squared displacement 
as well as other self-similar properties. In particular, the average motion remains 
invariant under the transformation {x,t) ^— {X^x,Xt), generalizing the self-similarity 
of the original Brownian motion, which is obviously recovered when if = 1/2 is set in 
Eq. ([2]). if is the self-similarity exponent of the process. When rewritten in terms of 



the Riemann-Liouville fractional integral operators (see Appendix A) fBm reads 

x(t)=a;o + oA'^''^'^%- (3) 

The corresponding propagator and diffusion equation of fBm have been derived in a 
number of ways in the literature (see, for example, [9] and [TO]). 

We will devote this paper to the Levy generalization of fractional Brownian motion, 
known as fractional Levy motion (fLm) |TT1 [121 [131 [Tlj . The Langevin equation defining 
the process is 

x(t)=xo + oA"''^'^""'a, (4) 

where ^a(t) is time-uncorrelated and distributed, for each t, according to a symmetric 
Levy distribution [15]. We recall here that symmetric Levy distributions are the 
symmetric solutions of the Generalized Central Limit Theorem and are parametrized by 
the stability index, a G (0,2], and the scale factor, a > 0. The characteristic function 
(i.e. the Fourier transform) of a symmetric Levy distribution La,cr{u) is 

^[L„,<,](A;) =exp(-a"|fcr). (5) 

In particular, for a G (0,2), the Levy distributions have algebraic tails, 

~ j^^^ W\ oo. (6) 



\U 



A Levy distribution with a = 2 is a Gaussian, 

1 r 



^^-^") = 2^^"H~4^J' 
and a is related to the second moment, {u^) = 2a^. 

The need for fLm originates in the observation of power-law (Levy) statistics for 
the displacements dx{t) in many physical systems of interest [ll[2], which is in contrast 
to the Gaussian character of fBm. In Eq. (Q, H is still the self-similarity exponent of 
the process. It can be shown [15| that finiteness requirements for certain moments of x 
restrict the admissible values of H to 
( (0,-] < a < 1 

Hel ^ ''^ - (8) 

[ (0,1] 1 < a < 2. 

Under these circumstances the s-th moment of x, with < s < a, behaves as 
{\x\') oct'^. 

Note that fBm is recovered from Eq. (Ill) if a = 2. Also, the increments dx{t) 
of the process are uncorrelated for H = 1/a. In this case, ordinary Levy motion is 
obtained, of which Brownian motion is a particular case {H = 1/a = 1/2). In the 
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present work we will compute in detail the propagator of fLm through path-integral 
techniques |T6], now familiar in both quantum field theory and statistical physics. 
Although the form of this propagator has been previously derived in a more abstract 
way by using self-similarity and stability arguments [11], the path- integral calculation 
offers a new insight which might help extend the range of applications of fLm and even 
tackle more comphcated situations. Two of us recently computed the propagator of 
fBm by path integral methods [T7| taking advantage of the fact that the measure can 
be written as the exponential of an action quadratic in the fields. As we will see the 
fLm path-integral measure is not Gaussian in the fields (except for a = 2) and the 
computation becomes quite more involved. 

The rest of the paper is organized as follows. In Section [2] we construct the 
appropriate probability measure on the space of realizations of the noise for fLm. The 
propagator is defined as a particular expectation value. Section [3] gives a detailed 
calculation of the propagator of fLm in the path integral formalism. In Section H] the 
fractional diffusion equation satisfied by the propagator of fLm is worked out. Section 
[5] contains the conclusions. Appendix A collects some basic definitions on fractional 
integrals and derivatives. 



2. Construction of the path integral measure 

Assume that the motion of a particle is defined by a stochastic differential equation. 
The propagator G{xt, T\xq, 0) is, by definition, the probabihty to find the particle at 
X = xt at time t = T ii initially, t = 0, it was located at x = xq. As mentioned above, 
the main objective of this paper is to compute the propagator associated to Eq. (jl]) (we 
drop the subscript a of ^ from now on) by means of path integrals. Consider trajectories 
x{t) : [0,T] M with boundary conditions x(0) = Xq and x{T) = xt- From Eq. (jll) 
we immediately deduce that the boundary conditions of x{t) are translated into the 
following constraint on ^ (t) : 

^D-''^'/-^'^ = XT-Xo. (9) 

The essential object in the path integral formalism is the probability measure 
V{^{t))V^{t) on the space of realizations of the noise, i.e. on the space of maps 
C,(t) : [0,T] M. Once it is constructed, the propagator is defined as the following 
expectation value: 

G(xr, T|xo, 0) = y 5 [oD^''^'^''-'^ - (xt - Xq)) Pmmit)- (10) 

In order to construct the measure associated to the Langevin equation (jlj) we 
will discretize the time in + 1 points tn := ne, n = 0,1,..., A^, with e := T/N. 
The continuum limit, N ^ oo, will be taken eventually. Each path is discretized 
according to x„ := x(t„). The appropriate discretization of the noise is made by 
taking ^(t„) = e^^+^/"^„, where each ^„ is an independent random variable distributed 
according to a symmetric Levy distribution of index a. The factor g^i+i/" ensures the 
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correct time-dependence of the finite moments of x, (|a;|^) oc t*^, < s < a. Therefore, 
the probabihty measure is naturally defined as 

N 

V{m)mt) = X{La.Ain)^in. (11) 

n=l 



Using the definition of the fractional integral, Eq. (lA.ll) . the constraint (19]) can be 
written as 

], [\t - rf~'/-ar)dr = Xr - x^. (12) 

T[H - l/a + 1) Jq 

Discretizing as prescribed above: 
/ (T - rf-y-ar)dT = ^ ^ne) / (T - rf-'^^dr (13) 

Jo J{n-l)e 



' n=l 

and we can easily solve for ^at in terms of n = 1, . . . , — 1: 

N-l 

^N = A-Y,BnU (14) 

n=l 

with 

V{H-l/a + 2). 
A '-{xT-xo), 

S„ := (Ar-n + l)^-i/"+^ - (Ar-n)^-i/"+\ n = 1, . . . , AT - 1. (15) 

Now, we are ready to write an explicit expression for the expectation value defining 
the propagator ffTOl) . Namely, 

G{XT, T|xo, 0) = ^iin /(T, N) f sUn-A+J^ Bn^n) J] L.A^dU, (16) 

\ n=l J n=l 

where f{T,N) is a normalization factor which will be determined at the end of the 
calculation. Equivalently, integrating over in Eq- (fT6|) : 

„ / N-l \ N-l 

G{XT, T\xo, 0) = Jhn /(T, N) / L,,. U " E H LaAU^^n- (17) 

\ n=l / n=l 

It is instructive to show that the path integral of ordinary Brownian motion, usually 
introduced in a different fashion, coincides with Eq. ( |T6l) when H = l/a = 1/2. The 
Langevin equation for Brownian motion is (recall Eq. ([1])) 

x{t)=Ut)- (18) 
The propagator is customarily introduced as 

1 '■^ 



G{xt, T|xo, 0) = / 5{x{0) - xo)S{x{T) - xt) exp ( -— J x{tydt ] Vx{t), (19) 
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where the paths are weighted by the classical action of the free particle. Now, one can 
choose the velocity, v{t) = x{t), as the integration variable. The transformation is linear 
and the Jacobian does not depend on the fields. The boundary conditions are simply 
translated into 

r-T 

V{t)dt = XT - Xq. (20) 





Therefore, we can write 

\2, 



v{tydt Vv{t). (21) 



< dt;,, (22) 



40-2 



G{xt, T|xo, 0) = y ^{^j^ v{t)dt - {xt - xo)^ exp (^-^ 
If we discretize the paths as above we get 

G{xt, T\xo, 0) = Jini^ /(T, A^) / 5 ( e ^ w„ - (xr - Xq) j JJ exp J^- 

\ n=l / n=l 

where /(T, N) is a normalization factor. Finally, with a last change of variables, 
'■= (■^^'^Vn (redefine f{T,N) as needed): 

G{XT, T|xo, 0) = jhn /(T, iV) /" 5 ( ^ - ) ([ exp (-^^n) d^, (23) 

\n=l ^ / n=l ^ ^ 

which is exactly Eq. ([T6D for = 1/a = 1/2 (recall Eq. 

3. Path integral computation of the propagator of fractional Levy motion 

The computation of f|T7j) is performed by repeated use of the identity 

L^A^)L^Ay - Ax)dx = (1 + A")-i/"L,,. ( (1 + A")V» ) ' ^ ^' (^4) 

where La^a{x) is a symmetric Levy distribution with index a and scale factor a, and 
A is a positive real number. The proof is straightforward. Define La^„{x) := La^„{\x). 
Then, 

/oo 
i^a,.(a;)L,,.(l/ - Aa;)dx = T-\T[L^,„]J'[L^,„]]{\-^y). (25) 
-oo 

Using that T[La^cr]{k) = \X\^^J^[La,a]{k/ ^) and JF[La,o-](/i;) = exp{—a"\k\") we get: 

^[L„,,](A;)^[I,,,](A;) = |A|-i exp(-(T"|(l + |A|-")i/"A:r). (26) 



And Eq. (|24|) follows easily. 

Let us go back to Eq. (IT7|) . Using Eq. (12^ we integrate out ^i: 
GixT, T\xo, 0) = Jhn /(T, AT) j L.,, (i+g)!/^ j H ^"><^(e")den. (27) 



Fractional Levy motion through path integrals 
Integration of ^2 yields 

/ 



G(xT,T|xo,0) = lim /(T,iV) / 



\ 



{{1 + Bfy/- (1 + ) 



l/a 



N-1 



n=3 



N-1 



= lim f{T,N) / 
And after — 1 integrations: 



(1 + 5f + B^y/^ 



71=3 



G{xT,T\xo,0) = lim f{T,N)L, 

N—*oo 



A 



(28) 



(29) 



It remains to compute 



lim 



A 



JV-l 



lim T{H-l/a + 2 f^„f'' x 



l + Y,{iN-n + l)^-V«+i _ (AT _ n)^-i/"+i)' 



n=l 



-l/a 



(30) 



In the sequel, g{N) ~ /'-(A^) will mean that g{N)/h{N) —>■ 1 when A^ — > 00. First 
observe that 

N-l 



n=l 



TV-l 



(i/ - l/a + l)°Ar"^'i ^ (1 



n=l 



iV 



where we have used 



n 1 
1 \ 

N N 



H-l/a+l 



1 

N 



n X H-i/a+i H -1/a + l f n\ 



N 



1 

A^ 



Now, we note that 
y-v n \ "^-1 1 

n=l 



NJ N 



1-1/N 

(1-m)°^-Mm 

l/N 



aH 



1 

A^ 



aH 



JSfaH 



and we have almost reached our goal. Combining Eqs. (131]) and (!33ll we get 

y ((AT - n + l)^-V«+i _ (AT - nf-^'^^'Y ~ (ff - l/a + 1)" ^,^^ 



Inserting this in Eq. (!30|) : 



lim 



a 



i7)i/"r(i^ - l/a + 1) 



Hence, 



G(XT, T|xo, 0) = /(T)L,,. ( {aHy/''T{H - l/a + 1) 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 
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and f(T) can be determined by normalization, G{xT,T\xo,0)dxT = 1: 

f(T) ^ i^^y-nH-l/a^l) ^ ,3,) 
SO that the final expression of the propagator is 

G[xt, T\xo, 0) = ^ L^^„ I [aH) ' T{H - 1/a + 1) 1 . (38) 

Summarizing, we have obtained that the propagator of fLm is a Levy distribution 
depending on the combination x/t^, so that the average motion is self-similar with 
exponent H. 

4. Fractional diffusion equation 

For the sake of completeness we derive in this section the fractional diffusion equation 
which governs the time evolution of the propagator of fLm. Denote by G{k,t) the 
Fourier transform of G{x,t\xo,0) with respect to x. Using Eq. ( !38l) . the form of the 
characteristic function of a Levy distribution, and the properties of the Fourier transform 
under rescaling: 

= ^''p(-(57rH&^'*'°)- '''' 

Differentiating with respect to t: 
Fourier inverting, recalling the definition Eq. flA.3p and the identity 



we find: 



d\x\ 



(k) = -\krf{k), (41) 



d a'^t"^'^ 8°' 

Therefore, the propagator of fLm satisfies a space-fractional diffusion equation with 
time-dependent diffusivity. Eq. (H2l) was recently derived by different methods in [T8] . 

The equation for the propagator of fBm (appeared in [10]) is obtained from Eq. (j^2l) 
in the particular case a = 2: 

which is a diffusion equation with time- dependent diffusivity. 

Finally, if H = 1/a = 1/2 we retrieve the standard diffusion equation associated to 
ordinary Brownian motion: 

^G{x,t) = a'-^G{x,t). (44) 
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5. Conclusions 

The Langevin equation defining iLm consists of two main ingredients: a time- 
uncorrelated stocliastic noise distributed according to a Levy distribution and a 
fractional integral operator which generates the time correlations. In this paper we 
have derived the propagator of fLm (which was deduced in [TT] by using self-similarity 
arguments) through path integral techniques. That is, we have explicitly constructed 
a probability measure on the set of realizations of the noise and precisely defined the 
propagator as an average over this measure space. The computation of the propagator 
has been performed by discretizing the paths and carefully taking the continuum limit 
at the end. The fractional diffusion equation associated to fLm has also been derived. 
We hope that the heuristic power of the path integral formalism will provide new insight 
on the calculation and help address more complicated cases. 
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Appendix A. Riemann-Liouville fractional integral and differential 



The books [191 EO] are excellent introductory texts to fractional calculus containing, in 
particular, the following definitions. 

Assume that / : M ^ M is a sufficiently well-behaved function. The Riemann- 
Liouville fractional integral operators of order a are defined as 



As for the Riemann-Liouville fractional differential operators of order a, the 
definition is 



where m is the integer number verifying m — 1 < a < m. 

Finally, the Riesz fractional differential operator is defined as the symmetric 
combination 



operators 






(A.2) 




a' 



-1 



(A.3) 
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